A VARIATIONAL METHOD FOR MULTISTAGE LAUNCH 


VEHICLE OPTIMIZATION 


by Fred Teren and Omer F. Spurlock 

Lewis Research Center 
Cleveland, Ohio 


TECHNICAL PAPER proposed for presentation at 


Third Aerospace Sciences Meeting of the American 
Institute of Aeronautics and Astronautics 
New York, New York, January 24-26, 1966 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 




A VARIATIONAL METHOD FOR MULTISTAGE LAUNCH VEHICLE OPTIMIZATION 


by Fred Teren and Omer F. Spurlock 
Levis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 


~o 

H 

i.O 


ABSTRACT 


\ATtO 


The methods of the calculus of variations are used to maximize payload 
capability for multistage launch vehicles. The method of solution uses the 
Lagrange multipliers to determine the optimum thrust direction profile, as 


well as to construct partial derivatives of payload with respect to the stage 
propellant loadings and a booster steering parameter. These derivatives are 
used to terminate the stages and/or as terminal equations to be satisfied. 
Maximum payload can thus be achieved with a single solution, rather than with 
a family of parametric results. Constant thrust and specific-impulse opera- 
tion is assumed for each upper stage (booster thrust and specific impulse 
vary with atmospheric pressure), and structure weight can be either fixed or 
a linear function of the stage propellant loading. Two-dimensional flight 
in a central , inverse- square gravitational field is assumed. 

Numerical results are presented for two- and three -stage launch vehicles 
flown to circular orbit and Earth escape, respectively. Parametric results 
are presented and compared with the overall optimum solution obtained by use 
of the variational technique. 

INTRODUCTION 

A problem that frequently arises in trajectory optimization studies is 
that of determining the maximum payload capability of a multistage launch 
vehicle flown to a prescribed set of burnout conditions. If all vehicle pa- 
rameters are specified, the problem reduces to that of finding the optimum 
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steering profile. In many cases, however, not all these parameters are speci- 
fied, and those left unspecified can be varied to maximize payload. 

A typical situation that occurs in the design of future launch vehicles 
is one in which the stage thrust levels and propellant flow rates (and, for 
practicality, the gross launch weight) are specified, but some or all the 
stage propellant loadings are left unspecified. The unspecified propellant 
loadings generally can be varied to achieve maximum payload capability for 
the vehicle. 

An additional optimizing parameter frequently is available in the booster 
steering program. Since the booster stage operates in the atmosphere , the 
booster thrust direction profile is shaped to minimize aerodynamic heating 
and loads and is not available for complete optimization. A single degree 
of freedom remains, however, corresponding to the magnitude of a short pitch- 
over phase following the initial vertical rise. This degree of freedom, some- 
times called the booster kick angle, determines the amount of trajectory loft- 
ing during boost phase. Since the upper stages operate essentially under 
vacuum conditions, the steering program for these stages is available for 
complete optimization. 

Many authors (e.g. , refs. 1 to 7) have treated the problem of optimizing 
the stage propellant loadings of multistage vehicles. None of these authors, 
however, has attempted to optimize the steering program for these vehicles. 
Others (e.g., refs. 8 to 10) have used the calculus of variations to optimize 
the steering program for various rocket vehicles. In particular, reference 11 
treats the problem of optimizing the steering program of a multistage launch 
vehicle. Reference 11, however, does not consider the problem of optimizing 
the stage propellant loadings or booster kick angle. 

Recently, Mason, Dickerson, and Smith (ref. 12) have considered the prob- 
lem of simultaneously optimizing the steering program and the stage propellant 
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loadings of a multistage launch vehicle. These authors followed the approach 
of Deribov (ref. 13) and Hunt and Andrus (ref. 14) in formulating the variational 
problem. 

The present report was written concurrently with reference 12 and presents 
a method which allows the propellant loadings, booster kick angle, and upper- 
stage steering program to be simultaneously optimized. The variational approach 
is somewhat different from that used in reference 12. By following the method 
of reference 15, the maximizing functional is written as the sum of the final 
payload and a constraint integral for each of the upper stages. The resulting 
boundary equations supply partial derivatives of payload with respect to the 
unspecified parameters. These derivatives are then used, along with the re- 
quired burnout conditions, as terminal equations to be satisfied. The analysis 
does not require that all the stage propellant loadings (or booster kick angle) 
be optimized. Equations are developed for optimizing payload with respect to 
any combination of unspecified parameters. 

The variational equations for optimizing vehicle parameters have been in- 
corporated into a digital computer program used previously at Lewis for para- 
metric launch -vehicle studies. Some of the procedures used to obtain numeri- 
cal results with this program are discussed in reference 16. Numerical results 
are presented for two- and three -stage launch vehicles flown to circular orbit 
and Earth escape, respectively, to demonstrate the feasibility of the varia- 
tional approach. Parametric results are presented showing the variation of 
payload with propellant loadings and booster kick angle. The resulting pay- 
load envelopes are then compared with the overall optimum points generated 
directly by means of the variational technique in order to verify the equations. 

ANALYSIS 

The problem to be solved is to determine the maximum payload capability 
of an N-stage launch vehicle flown to a specified set of burnout conditions. 
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The analysis admits atmospheric effects during booster phase but assumes vacuum 
operation for all other stages. Because of these atmospheric effects, the 
booster steering program is assumed completely specified (e.g., zero angle of 
attack), except for the booster kick angle. The upper-stage steering program, 
however, is unconstrained and is determined to maximize payload. The calculus 
of variations is used for this purpose. 

Each of the upper stages is assumed to operate at a fixed (constant) value 
of thrust and propellant flow rate. These values may be zero, so that coast 
phases are admitted. The structure mass for each stage is assumed to be a 
linear function of the stage propellant loading, defined by 

m s = m R + kmp 

where m s is the total structure mass, m H is the fixed mass, m p is the stage 
propellant mass, and k is the propellant sensitive mass fraction. (All sym- 
bols are defined in appendix A. ) 

In addition to the variational trajectory, provision is also made for add- 
ing an additional velocity increment AVj after the desired orbit conditions 
are achieved. This velocity increment is achieved by use of the final stage 
for propulsion. ‘The amount of propellant required for this maneuver is calcu- 
lated by use of the standard impulsive velocity equations. 

Variational Problem 

Since the booster steering program is not subject to complete optimization, 
the booster stage is not treated in the following Euler -Lagrange equations. 

The booster degrees of freedom (propellant loading and kick angle) are included 
by allowing variations in the position and velocity at second-stage ignition. 

The associated equations, along with the equations for optimizing upper-stage 
propellant loadings, are treated in the boundary equations resulting from the 
variational analysis. 

The variational problem to be solved is that of finding the upper -stage 
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thrust program which maximizes the payload capability of an N-stage launch 
vehicle for given boundary conditions. This problem can be formulated as a 
generalized Bolza problem. By following the treatment in reference 15, the 
functional to be minimized can be written as 


T ~ 


= g + 



i=2 J ti_i 


F ± dt 


( 1 ) 


where g is a function of initial and final conditions to minimized and 
consists of a set of constraints to be applied to each of the upper stages, 
added together with the aid of Lagrange multipliers. For this problem, the 
payload is to be maximized, so that 


S = " m PL 

The constraint equations are 

9 Ti 

f - u + J± — orr - — sin \lr = 0 

11 _2 m 


f o • = rd) + 2uo) - — cos = 0 
2i m T 


f 3i = r - u = 0 
f 4i =q)-^ = 0 
^5i ~ e. + 3i = 0 
applicable on the interval 


ti_i £ t £ ti 


( 2 ) 

(3a) 

(3b) 

(3c) 
(3d) 
( 3e ) 


for i ~ 2 , 3 , . . . , N 

Equations (3a) to (3d) are the two-body equations of motion written in two- 
dimensional polar coordinates with an inverse square force field acting. The 
thrust direction \\r and the state variables are defined in figure 1. Equa- 
tion (3e) defines the propellant flow rate. 

Equations (3) are combined to give 


E. = 


- 22 Li f ji 


i = 2, . . . , N 


(4) 


5 



where Ajj_ are undetermined Lagrange multipliers, which are functions of time 
since the constraint equations must be satisfied at all points of the trajectory. 

Euler -Lagrange Equations 

As shown in reference 17, a necessary condition for g to be minimized 
is that the Euler -Lagrange equations be satisfied. The Euler -Lagrange equa- 
tions are 



where x. are the problem variables 
J 

x 1 (t ) = u > 

X 2 (t) = 03 

x 3 (t) = r 

, + , f (6) 

X 4 (t) = cp 
xs(t) = m 
x g(t) = \| i j 

The Euler -Lagrange equations for the present problem can be written explicitly 
by use of equations (3) and (4): 


~ ScuAg — A^ 

(7a) 

u ^4 

= -2aiA-i + — Ap - — 
1 r * r 

(7b) 

= - + cb^A^ + cbAg 

(7c) 

= 0 

( 7d) 

T 

= — o (A-i sin \|r + Ap cos 

(7e) 

T 

cos \|r - Ag sin = 0 

(7f) 

s (7) (and. subsequent equations) 

apply separately to each of the upper 


stages. The subscript i has been omitted for simplicity. 

Equation ( 7f ) determines the thrust direction (for T ^ 0): 
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(8a) 


^1 

tan i\i = — 

A 2 

sin \| i = — - --- -- ■ (8b) 

aAi + 4 

A 2 

cos \J| = — , . (8c) 

The choice of sign in equations (8b) and (8c) has been determined by use of 
the Weierstrass E-test (ref. 18). 

Equations (7a) to (7d) must be integrated, along with the equations of 
motion to determine the thrust direction and the optimum trajectory. It is 
shown later that equation (7e) need not be integrated. Equation (7d) is 
easily integrated to give 

= Constant (9) 

First Integral 

Since the function F does not contain the independent variable (time) 
explicitly, a first integral to the Euler -Lagrange equations exists (for each 
stage), which can be stated as (ref. 17) 



( 10 ) 


Equation (10) can be written explicitly by use of equations (3) and (4): 

C - - a?rj Aq - 2 ucdA 2 + u ^3 + ^4 " + “ (Aq sin ^ + A 2 cos = 0 

(ID 

Equation (11) is used later in the analysis to eliminate both C and A5 
from the boundary equations. 

Weierstrass -Erdmann Corner Condition 
The boundary conditions on the Lagrange multipliers at staging points 

can be derived with the aid of the Weierstrass -Erdmann corner condition 

c)F 

(ref. 17). This condition states that (j = 1, . . .,6) must be continuous 
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at such corners. For the present problem, this condition implies that all the 
multipliers are continuous at staging points, hence continuous throughout the 
trajectory. 


Trans versality Equation 

The relation between changes in boundary conditions and changes in J is 
expressed by the general transversality equation (ref. 17). For this problem, 
the transversality equation can be written 



This equation can be written explicitly by using the definition of F and the 
first integral 


1 + dg (13) 

i-1 

The subscript i has been used with C and dm since these variables may 
be discontinuous at staging points. 

Boundary Equations 

If some of the problem variables (state conditions or control variables) 
are not specified, values should be chosen which minimize J (or, equivalently, 
maximize payload). 

According to reference 17, minimizing J is accomplished by setting dJ 
equal to zero. Equation (13) has the form 
m 

Gj dxj = ° ( 14 ) 

If the m problem variables x, are all independent, dJ wil vanish if, and 

only if, each term on the right side of equation (14) is independently set equal 

to zero. For specified variables x., the allowable variation dx- is zero: 

J J 9 

for unspecified Xj, the coefficient Gj must loe set equal to zero. Equa- 



dJ = 



i= 


(c ± dt + A-j_ du + rAg doo + A3 dr + A^ dcp + A5 


8 



tion (14) can thus he Interpreted as a total differential of J , and 


r aj 


Equation (13) is not suitable for this interpretation, since the variables are 
not all independent. In the following section the dependent variables are 
eliminated by expressing the dependence explicitly. 

Consider first the terms in equation (13) involving the variations of the 
state variables. 


(a-^ du + rAg dm + A 3 dr + A 4 dp).^ ^ £(A-^ ^ u ) + “ (A]_ du) + (rAg dm) 

- (rAg dm) + (A 2 dr) + - (A 3 dr) + (A 4 dp) - (A 4 dp) 1 


+ (A]_ du + rAg dm + A 3 dr + A 4 dp) - (A x du + rAg dm + A 3 dr + A 4 dp) 

t = t - 


'N 


t=t 


(15) 


where the superscripts - and + refer to conditions before and after staging, 
respectively. Since the state variables and Lagrange multipliers are continuous 
throughout the trajectory, the summation on the right side of equation (15) is 
identically zero. At t = t-^, the variations in state conditions are due to 
the allowable variations in booster burning time and kick angle. 


(dxj 

J t=tp 


dx. dx* 

r— 1 “• dT-i + sr— 1 - da 
dx n x aa 


j = 1, • 


(16) 


where a is the booster kick angle and is the booster burning time. (in 

general, is the burning time of stage i, that is, T j_ = ^or 

generality, the state variations at t = t^ are expressed in terms of a set 
of generalized (independent) state variables, r^, k = 1, . . . , 4, so that 


dx-; 


E 

k=l 


Sx. 




j • i 


) 


4 


(17) 


Combining equations (15) to (17) yields 
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i — * + 



^ (\ du+rA 2 doj + A 3 dr + A 4 <%)** ± = ^ 

1_ 1_ k=l ' k 


+ rA 2 §^ +A 3 |t +A 4 ^ 


d \/ 


‘ + rA 2 5f + A 3 s| + A 4 




t=t» 


d.T i 


(18) 


Inasmuch as a suitable form for the state variation terms in equation (13) 
has been obtained, the mass and payload variation terms are now considered. 


Since the propellant flow rate is constant for each stage, the variables m 4 
can be expressed as functions of the burning times of the stages. Specifically, 


4 

-*-s 

j^ 1 + 

+ m H,2j 

A 


(19a) 

4 

II 

3 

i 

j(l + k 2 )p z x z 

+ m H ,Z 

“ Pi T i 

J 

> i = 1, . . . , N 

(19b) 


where m Q is the lift-off mass. The summations in equations (19) (and all 
other similar summations) are defined to be zero whenever the lower summation 
limit exceeds the upper limit. The superscripts 0 and f refer to condi- 
tions at the beginning and the end of each stage, respectively. 

In order to calculate the final payload, the velocity impulse discussed 
earlier must first be considered. This velocity impulse is added after the 
desired orbit is achieved, with the use of the N th stage for propulsion. After 
this maneuver is completed, the stage hardware is jettisoned. The propel- 
lant required for the velocity impulse is 


Am m f (-\ “ AV I P N/%\ 

^1 = - e ) 

By use of equation (19b), the final payload can he written 


m. 


PL 


T {^° " 0 1 + + m H,z| ' |3 N t n| " 


k N% T N m H,N 


where y has been defined as 

r - -k„ + (1 + k N )e- AV I P » /T n 


( 20 ) 


( 21 ) 



can now be expressed in terms of the 


0 f 

The variations diriY, dmi , and dirip^ 

variations in burning times by differentiating equations (19) and (20), 

(1 + k z )|3j dTj | (22a) 


dm9 = - 






1=1 

w-i 


(! + k z )p z dT Z - pi dTi 


X — 1 ^ ^ IT 


dffipL = " r ^ ( 1 + k pPj dx l - (r + %)P5r dx N 


(22b) 


(22c) 


By use of equations (22), the terms in equation (13) involving dn^ can be 
simplified and expressed in terras of the variations dr. 


N-l Ebl 

(A 5 dmi^i^ = g ^5 k iPi dT i ’ A 5 ^ + k pPi dT i 


*Vi 

-AgPjj dTjj + A 5 (l + k f) Pq dx^ 


(23) 


where 


As = (Ac) 


5 t=ti 

Also, the time variation terms in equation (13) can be expressed in terms of 
variations in the stage burning times: 


N 


N 

jQ (q dt)^_ x - g C t dTi 


(24) 


The variation of g is calculated by combining equations (2) and (22c): 


1=1 


dg = r E (1 + k i )p i dti + dr ; 


N 


(25) 


The variable A^ can be eliminated from equation (23) by use of the following 
identity: 


$ ♦ £(£- > i - 1 ) + it <*t - ^ * * 

7.=i+i vl / TZ T+l V l ' i=i+l 


R+l 

M° 


( 26 ) 




p z =° 


It is convenient to define 
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(27a) 


S i = — ' A 5 
1 Pi 5 


Q 0 i ,i-l 

b_- = - — - Act 

1 Pi 5 


sf = o 


s? = o 


Pi i o 




i = 2, . . ., N 


-.f 




Pi t o 


Pi = o 


Pi = 0 


s " = " ft + rA 2 5^; + a 3 5% + a 4 


tj 


Also, note that Ag = 0 for coast phases, so that for It 

1 _ i-1 
A 5 - A 5 

By use of equations (27) and (28), equation (26) becomes 


= 0, 


»S- 



(sf - 3°) + >| 


(27b) 

(27c) 

(27d) 

(27e) 

(28) 

(29) 


Equations (16) to (18), (23) to (25) and (29) are now combined with equa- 
tion (31) to give 
N-l 


[ (1 + k i )fJ i[ r + ,2^ (s i • s °] + Pi s ij a n 

eI?o 


r— t j 

C i dT i + [(v + %)Pn + PN^| dT N “ ( 

Pi=0 


h if + rAz s + A3 s + A4 sf) t=ti da 


3=1 V J 


+ rAc. ^ r ~ + A* + A- &L.\ 


5o) 

x 2 Stj, 


3 T " 4 a ij l " ,J 


d Tll = 0 


(30) 


t=t 


N 


The variations in equation (30) are all independent, so that the form of equa- 
tion (14) has been achieved, with 


G(xi) = (1 + k i )p i Jy + - S°)j + PiS?, p ± 


t 0 

i = 1, . . . , W - 1 


(31a) 
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G(t ± ) = C ± , pi = 0 1 = 2, 

G(tjj) = (t + + 3jjS^ 


, N - 1 




1 S + ^ A 2 S + A 3 s 


, , cir .. 
+ A* x~ + A 


A 

x /t=t x 

G(, 'o 1 ‘ (\ 5^7 + rA 2 + *3 3^ + X 4 t^) 

t — t 


J = 1, 


(31b) 

(31c) 

(31d) 
, 4 ( 31e ) 


From equations (ll) and (27 ) , 

fe-*) 


T, 

A-, + 2uodA ? - uA., - (bA. - 

x ^04 m i 


r£ 

D i 




i + A S 


J t=t i 


Pi 


S? 


C i = 


2 + 2uajAg - uA^ - aoA^ 

m i 


- tv?: 
r 


-%4 


, 9 9 

Ai + *2 


”* t='b i _ 1 


Pi 


K -% “ a/or ^A-i + 2uo)Ao - u?w - cjdAJ 
r Z / 1 ll p 1 =0;t i _ 1 <t<t i 


(32a) 

i = 2, . . . , N; 

P ± / o 

(32b) 

(32c) 


Since equations. (32) do not contain C or A^, equation ( 7e ) need not be inte- 
grated to evaluate equations (3l), as indicated earlier. 

Boundary Value Problem 

The determination of an optimum trajectory requires the simultaneous 
integration of the equations of motion and the Euler -Lagrange equations. A set 
of initial conditions (state variables and Lagrange multipliers) and staging 
times is required in order to specify the trajectory uniquely. 

The trajectory thus generated must satisfy N + 5 final conditions, corre- 
sponding to the N + 5 independent problem variables in equation (30). For 
specified variables, the final conditions have the form 


x j 


(33) 


where the subscript d indicates the desired final value. For unspecified vari- 
ables, equations (31) supply auxiliary final conditions with the form 
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G(xj) =0 (34) 

Some of equations (33) are easily satisfied; for example , a specified "burn- 


ing time for any stage can be achieved simply by terminating that stage at the 
proper time during the integration. An iteration is required in order to satisfy 
the nontrivial final conditions , and variable initial conditions (equal to the 
number of final conditions) must be available. For the present problem, the 
Lagrange multipliers (A-p i = 1, . . .,4), the burning times of all stages 
being optimized, and the booster kick angle (if it is being optimized) are avail- 
able as variable initial conditions. 

The size of the iteration loop can be reduced by using the fact that equa- 
tions (7) are homogeneous in the A T s. This implies that the choice of any one 
A is arbitrary and serves only as a scale factor for the others. The value of 
this multiplier can be chosen to satisfy one of equations (31). 

The iteration size can be further reduced when A^ = 0 is a required final 
condition, which occurs when the travel angle is unspecified (ref. 8). Since A 4 
is a constant and is continuous across staging, this final condition can be satis- 
fied at t = t-j_; A^ is thus removed from the iteration. Equation (31d) can also 
be evaluated at t = t-^, and the booster kick angle optimization can be similarly 
removed from the iteration. 

The burning times of optimized stages can sometimes be removed from the 
iteration, along with an equal number of equations (31) ♦ The principle involved 
is similar to that used in removing equation (31d): If all variables in the 

equation G(t ^) =0 can be calculated at (or previous to) t = t^, stage l can 
be terminated (during stage l) whenever the optimizing equation is satisfied. 

The number of optimized stages which can be thus removed depends on the charac- 
teristics of the particular problem. 

Consider first two powered stages and r m , l < m, with k^ f 0. Equa- 

tion (31a) must be set equal to zero for i = l and i = m, and the two result- 
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.rig equations are combined to give 
m-1 

TTT j + (sf ' s?) + r 


qf 

+ V ^ 


sS = o 


i=2+l 


(35) 


Equation (31a) (for i = 2) can be satisfied by the choice of one of the Vs. 
Equation (35) contains terms which can all be calculated at, or prior to, stage 
m cutoff, and this equation can be used to terminate stage m. Specifically, 
stage m is terminated when 


S 

m 


, f _ 1 + k m 


SP - 
m 


1 + ki 



(36) 


•P > . 

If k^ = 0, the term in disappears from equation (35), and this equation can 

then be used to terminate stage m-1: 




m-1 

m-2 

qf 

m-1 

q0 _ 

m 

t 

Z <-r 



i=2+l 

i=2 +1 


rtf 



S° - - 


- 0 

m = l + 1 


1 + k., 


m > 2 +1 


m 1 + k 


(37a) 


(37b) 


'The terms on the left side of equations (37) are evaluated at stage m-1 cut- 
off, and the terms on the right side are evaluated during previous stages. 

For the special case m = N, kjj ^ 0, equation (31c) must be set equal to 
zero. Equations (31a) (for i = l) and (31c) are combined to give 

,f N-l 


4 


y + k- 


•N 




r+T^ - $ + Z ( s ( - s °i) 


i=2+l 


(38) 


Equation (38) is used to terminate stage N. 

For coast phases to be optimized, equations (31) supply criteria for termi- 
nation of the previous stage, similar to equation (37). For such cases, stage 
7-1 is terminated when C 7 = 0, where stage 2 is the coasting stage to be 
optimized. 
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RESULTS 


In order to demonstrate the validity of the equations and the feasibility 
of the variational technique , parametric results are presented and compared 
vith the overall optimum solutions obtained by using the variational technique. 
Two- and three-stage launch vehicles were optimized by use of the equations 
developed. The results presented include a two-stage vehicle flown to circular 
orbit and a three-stage vehicle flown to Earth escape through a circular parking 
orbit. In these results, both fixed and variable hardware weights were used, 
and the propellant loadings and booster kick angle have been optimized. Para- 
metric results are also presented, which show the variation of payload with pro- 
pellant loadings and kick angle. 

Vehicle Definition 

The vehicle chosen for this study is a hypothetical three-stage launch 
vehicle consisting of two chemical stages and one nuclear stage. The assump- 
tions on propulsion and weights are listed in table I. 


TABLE I. - LAUNCH VEHICLE PROPULSION AND WEIGHT DATA 



Stage 



First 

Second 

Third 

Thrust , lb 

7.5X10 6 (Sea level) 

1.5X10 6 

2.5X10 5 

Specific 

264 (Sea level) 

428 

850 

impulse, sec 
Fixed 

305 (Vacuum) 
245,000 

70,000 

35,000 

hardware 
weight, lb 
Propellant 

0.030 

0.033 

0.120 

sensitive 
fraction 
Drag reference 
area, sq ft 

855 




The first stage is based on the Saturn SI-C stage, consisting of five F-l 
engines using RP-LOX propellants. The engine performance data are used for 
illustrative purposes only and are not necessarily consistent with present 
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31-0 values. The second stage consists of one M-l engine using liquid hydrogen 
and liquid oxygen as propellants. The third stage is a nuclear stage. The 
launch thrust -to -weight ratio was fixed at 1.25 for this study , and the launch 
azimuth was 90°. 

Numerical Results 

A typical propellant tank sizing study was conducted with the vehicle de- 
fined in table I. In this study, the stages are assumed to have variable tank 
size, and the optimum propellant capacities are determined by flying possible 
missions of interest. 

The first mission is flown to Earth escape energy with three stages. This 
mission is typical of lunar and planetary probes and orbiters. A parking orbit 
ascent mode is used, wherein all three stages are used to enter a parking orbit 
of a 121 nautical-mile altitude (similar to the Apollo mission), after which 
the third stage is burned to escape. The second burning of the third stage is 
assumed to be impulsive, with a velocity increment of 10 563 feet per second. 
The maximum payload capability for this mission is obtained by optimizing the 
booster kick angle as well as the propellant loadings of the three stages. By 
use of the equations derived earlier, maximum payload can be obtained with a 
single (converged) solution, represented by the optimum point in figure 2. 

With the parametric procedure, the maximum payload is obtained as the 
envelope of payload points with ail possible combinations of propellant load- 
ings and booster kick angle. The increased effort (number of solutions) re- 
uuired for this procedure is obvious. In addition to the parametric curves 
shown in fig'rre 2, each point on the curves had to he obtained at an optimum 
kick angle, which required an additional family of curves (not shown in the 
figure), with kick angle as the independent parameter. As can be seen from 
figure 2, the payload capability and optimized propellant loading obtained 
from the variational procedure agree with the parametric results. 
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Another mission of interest is two stages flown to a circu lar orbit, with 
the use of the first and second stages from table I. Results for this mission 
are presented in figures 3 and 4. In figure 3, with the use of the parametric 
procedure, payload capability is presented as a function of first-stage pro- 
pellant loading for various booster kick-angles. The maximum payload for each 
kick angle is determined from the figure and presented as a function of kick 
angle in figure 4. The case flown by use of the variational technique is repre- 
sented by the optimum point in figure 4, and the payload obtained is in agree- 
ment with the envelope, as before. 

Frequently, the propellant tanks of the various stages are sized for one 
mission and fixed at these values for all other missions. Propellant loadings 
can still be optimized in such cases but with the restriction that the propel- 
lant loadings cannot exceed the propellant capacities of the tanks. 

In figure 5, the three-stage Earth escape mission is reoptimized with 
fixed tanks for the first and second stages. The tank weights and propellant 
capacities used are based on optimum values for the two-stage orbit mission. 
Since the maximum propellant capacities of the first and second stages were not 
exceeded in the optimum case, the propellant loadings of these stages were off- 
loaded to the optimum values. The parametric results and the variational point 
are compared (as in fig. 2), and the results are in agreement. 

CONCLUDING REMARKS 

A technique is presented which allows simultaneous optimization of the 
thrust direction profile and vehicle control parameters for multistage launch 
vehicles. The agreement of parametric and optimum results presented in fig- 
ures 2 to 5 demonstrates the correctness of the optimizing eq ua tions. 

The amount of effort (and computer time) saved by the variational technique 
in determining the maximum payload capability can be seen by referring to the 
figures. In figure 2, for example, approximately 80 parametric data points are 
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required to optimize the three -stage propellant loadings and booster kick angle 
(not shown in the figure). When good initial guesses are available , the com- 
puter time required to obtain the overall optimum solution is about the same 
as for each parametric solution. The time saving is therefore proportional 
to the number of parametric cases required fctr complete optimization. 

APPENDIX A 
SYMBOLS 

C constant of integration 

E Weierstrass excess function 

e eccentricity 



f constraint equation 

G dJ fdx 

g function of initial and final conditions to be minimized 

J functional to be minimized by variational methods 

k propellant sensitive mass fraction 

m mass, slugs 

r radius , ft 

S functions defined in eqs. (27) 

T thrust , lb 

t time, sec 

u radial velocity, ft/sec 

v velocity, ft/sec 

x problem variable 

a booster kick angle , rad 

p mass flow rate, slug/sec 

X function defined in eq. (21) 

q generalized state variable 
A Lagrange multiplier 
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M- Earth force constant, cu ft /sec 2 

x burning time, sec 

9 polar angle, rad 

thrust direction, rad 
co angular velocity, rad/sec 

Subscripts : 
d desired value 

H fixed hardware 

I impuls ive 

i stage number 

j variable number 

k variable number 

l stage number 

N last stage 

EL payload 

p propellant 

s structure 

0 initial 
Superscripts : 

f end of stage 

1 refers to t = t-j. 

0 beginning of stage 

• derivative with respect to t im e 

+ after staging 

before staging 
-*■ vector 
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Figure 2. - Payload capability as function of second-stage propellant 
loading. Three stages to Earth escape via 121-nautical-mile park- 
ing orbit; variable tanks; first-stage propellant loading at optimum 
point, 4 028 000 pounds; third-stage propellant loading at opti- 
mum point, 317 000 pounds. 




First-stage propellant loading, lb 

Figure 3. - Payload capability as function of first-stage 
propellant loading. Two stages to 121-nautical-mile 
orbit. 
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Figure 4 - Payload capability as function of booster kick angle. Pro- 
pellant loadings optimized; two stages to 121-nautical-mile orbit; 
first-stage propellant loading at optimum point, 3 980 000 pounds; 
second-stage propellant loading at optimum point, 1 276 000 pounds. 



Figure 5. - Payload capability as function of second-stage propellant 
loading. Three stages to Earth escape via 121-nautical-mile park- 
ing orbit; fixed tanks; first-stage propellant loading at optimum 
point, 3 964 000 pounds; third-stage propellant loading at optimum 
point, 310 000 pounds. 
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